Dumb-bell swimmers 
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We investigate the way in which oscillating dumb-bells, a simple microscopic model of apolar 
swimmers, move at low Reynold's number. In accordance with Purcell's Scallop Theorem a single 
dumb-bell cannot swim because its stroke is reciprocal in time. However the motion of two or 
more dumb-bells, with mutual phase differences, is not time reversal invariant, and hence swimming 
is possible. We use analytical and numerical solutions of the Stokes equations to calculate the 
hydrodynamic interaction between two dumb-bell swimmers and to discuss their relative motion. 
The cooperative effect of interactions between swimmers is explored by considering first regular, 
and then random arrays of dumb-bells. We find that a square array acts as a micropump. The long 
time behaviour of suspensions of dumb-bells is investigated and compared to that of model polar 
swimmers. 

PACS numbers: 47.63.Gd, 47.63.mf 



INTRODUCTION 

Bacteria swim in a very different manner to macro- 
scopic animals because, at micron length scales, inertial 
effects are negligible compared to viscous forces. In this, 
the zero Reynolds number limit, the Navier-Stokes equa- 
tions reduce to the Stokes equations, which are time re- 
versal invariant. This leads immediately to the Scallop 
Theorem that any swimming stroke must be non- 
reciprocal - look different when run forwards or back- 
wards in time - if it is to result in a net motion. Examples 
of strategies developed by microorganisms to overcome 
the Scallop Theorem are rotating flagella, waving cilia 
and surface waves. Moreover initial experiments have 
shown that it is possible to design tiny swimming robots 
which may be developed for drug delivery or manipulat- 
ing payloads in microchannels [^. 

To gain insight into the generic behaviour of micro- 
scopic swimmers several model systems have been de- 
fined and studied both analytically and numerically. 
These include squirmers [3, |3| , Purcell's three link swim- 
mer 0, 0] , three linked spheres 0, 0] and the 'push- 
mepullyou' swimmer Q. All of these systems are po- 
lar swimmers, in that they undergo a non-reciprocal 
stroke, which picks out a preferred direction in time and 
hence results in swimming. Indeed, the Scallop Theo- 
rem suggests that apolar models should be immediately 
discounted as they will not swim. 

However, as pointed out by Koiller et al although 
this is true for a single apolar swimmer, a collection of 
such swimmers can swim. If their swimming strokes are 
not in phase, the motion taken as a whole is not recipro- 
cal in time, and motility is possible. The physical mech- 
anism which turns the motion of an apolar swimmer into 
a net displacement is the hydrodynamic forces between 
the swimmers. 

In this paper we study the simplest model of an apolar 



swimmer, oscillating dumb-bells. We give an analytic 
solution for the motion of a pair of dumb-bells, in the 
Oseen tensor limit, and valid for large separations. By 
numerically iterating these results we are able to gain 
insight into how the orbits of two swimmers depend on 
their relative positions, separations, and phases. We then 
consider regular arrays of swimmers and, in particular, 
suggest a way in which tethered swimmers could be used 
as a micropump. 

We next present results for suspensions where the 
dumb-bells are initially placed randomly with random 
relative phases. The collective behaviour of such groups 
of self-propelled organisms has received a lot of interest 
in the past decade. Remarkably, the dynamics of sys- 
tems as seemingly unrelated as flocks of birds, swarms 
of bacteria and vibrated rods can be treated within the 
common theoretical framework of active fluids 
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Within this class a distinction has been made between 
systems where the order is polar, for example, flocks of 
birds, and those with apolar order, such as vibrated rods. 

Since a dumb-bell is intrinsically apolar, a suspension 
of them provides a minimal microscopic model of an ac- 
tive apolar fluid. We present results comparing the be- 
haviour of a suspension of interacting dumb-bells to that 
of a suspension of three-linked-sphere swimmers 0, [1], 
which is a simple microscopic model of an active polar 
fluid. 



HYDRODYNAMIC INTERACTIONS BETWEEN 
OSCILLATING DUMB-BELLS 

We first describe the motility of oscillating dumb-bells 
at zero Reynolds number. Each dumb-bell comprises two 
spheres, of radius a, joined by a thin, rigid rod, whose 
length varies sinusoidally as D -|-^ sm{Lut). We stress the 
importance of the relative phase of the two dumb-bells: 
if they oscillate in phase or tt out of phase then their 
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combined movements remain reciprocal and the Scallop 
Theorem prevents any net motion However, for other 
values of the relative phase, motion will occur. To show 
this we use the Oseen tensor formulation of hydrodynam- 
ics, valid in the limit of zero Reynolds number [l3|. We 
present anaytic calculations, valid for large separations, 
and numerical results, valid for separations down to of 
order the size of a dumb-bell. 

Linearity of the Stokes equations which govern zero 
Reynolds number flows allows the fluid velocity to be 
written as a linear combination of the forces, f , acting on 
the fluid due to the motion of the swimmers 



u(x) = ^[G^(x)fA + G^(x)f 



(1) 



Here the subscript A labels the dumb-bells and the super- 
scripts 1,2 label the spheres of an individual dumb-bell. 
For spheres whose radii a is small compared to their sep- 
aration the Greens function, G^(x), may be taken to be 
the Oseen tensor. 
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The motion of the dumb-bells 



through the fluid is then determined by three ingredients; 
consistency of the fluid flow with the change in shape, i.e., 
the change in rod length, the constraint that each dumb- 
bell is force free, -I- = 0, and the constraint that 
each dumb-bell is torque free, (x^ — x^) Ai\ — 0. 

A short calculation then leads to an expression for the 
forces 
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-G|(xi) + G|(xi)}fij, 



(3) 



where terms of 0[{a/D)'^] have been omitted. The unit 
vector riA gives the direction of sphere 2 relative to sphere 
1 and thus describes the orientation of the dumb-bell. 
Applying the torque free constraint determines the an- 
gular velocity, f2^, from which we obtain an equation for 
the evolution of the dumb-bell's orientation 
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(5) 



Finally, the translational motion of each dumb-bell is 
given by 



dx^ 
dt 
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(6) 



where x^ denotes the 'centre' of the dumb-bell, x^ : — 
(xJi -I- x^)/2. As expected, the motion described by 
Eqns. ^ and ([6]) arises solely through interactions with 
other dumb-bells. 

We assume that the dumb-bells are in a dilute suspen- 
sion so that the separation, r, of any given dumb-bell 
from its nearest neighbour may be assumed to be large 
compared to its size, D. Under such circumstances the 
interactions between the dumb-bells may be expanded in 
a power series in {D/r) and only the leading contribu- 
tions retained. Integrating over a complete cycle leads to 
expressions for the changes in position and orientation of 
the dumb-bells after a single swimming stroke: 
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2nA 



{■o-A ■ Vba) + 2(nA ■ nB)(n_B • vba) 



5(nA • vba){t^b ■ iBAf | , 



where yba is the position vector of B relative to A. The 
amplitude A is given by 

A = nD^U^B sin(rysA) + sin(2?7BA) , (9) 

where rjBA is the phase of -B's swimming stroke relative to 
A's. The long time behaviour of a group of dumb-bells 
may be determined by numerically iterating Eqns. ^ 
and ([5]) to find the new positions and orientations of all 
the dumb-bells after each swimming cycle. 
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TWO DUMB-BELLS 

We first consider two dumb-bells lying in the yz— plane 
and both oriented along the z— direction. Dumb-bell A 
is initially at the origin and dumb-bell B is placed on a 
circle of radius AD centred on the origin. For two dumb- 
bells, varying the relative phase does not lead to any 
qualitative changes in behaviour and therefore we con- 
sider only the case r]BA = 7r/2. The sole free parameter 
is the angle, 9, that the position vector of B makes with 
the y— direction. As this angle is varied the hydrody- 
namic interactions between the dumb-bells change lead- 
ing to different long time behaviour, which we illustrate 
in Fig. m 

The predominant behaviour is for the two dumb-bells 
to adopt a perpendicular configuration in which one 
dumb-bell is oriented parallel and the other perpendicu- 
lar to their relative position vector. This stable arrange- 
ment appears at long times for all initial configurations in 
the regions labelled P in Fig.[T] An exemplary time series 
showing how the perpendicular configuration is reached 
is shown in Fig. [Hi). The stability of this state may be 
seen from a linear stability analysis of Eq. ([7]) . 

Similarly, linear stability analysis reveals that the rota- 
tional fixed point at = is unstable. However, small de- 
viations away from this fixed point do not lead smoothly 
to the stable perpendicular configuration. Instead there 
is an initial period during which the dumb-bells tumble, 
often several times, before they finally settle down, as il- 
lustrated in Fig. [T]^ii). Exactly at the fixed point, and in 
the absence of any fluctuations, the dumb-bells undergo a 
pure translational motion, swimming cooperatively in the 
direction of the dumb-bell with positive relative phase. 

The fully aligned configuration with 6 — tt/2 is also 
a rotational fixed point and again, exactly at this angle, 
the motion is purely translational with both dumb-bells 
moving in the same direction. However, in this case, the 
fixed point is a centre and small deviations away from it 
lead to an oscillatory cooperative motion, an exemplary 
time series of which is shown in Fig. [IJiii). This oscilla- 
tory motion occurs throughout the regions labelled O in 
Fig. [1] and is separated from the P regions by an addi- 
tional rotational fixed point at w 65°. 



FIG. 1; (Colour online) Orbits of a pair of coplanar dumb- 
bells. The two dumb-bells are initially parallel with A at 
the origin and B placed at the corresponding point in the 
diagram. Several regimes of long time behaviour are found. 
In region P the two dumb-bells adopt a stable perpendicular 
configuration as exemplified by (i) . This type of behaviour is 
also found in the shaded regions, however, before acquiring 
the perpendicular configuration the two dumb-bells tumble 
for a period of time, as shown in (ii). In the regions labelled 
O the dumb-bells move along parallel trajectories with an 
oscillatory motion, as illustrated in (iii). Finally, along each 
of the dashed lines there is no rotational interaction and the 
dumb-bells undergo a pure translation along straight parallel 
paths. 



LINEAR CHAINS OF DUMB-BELLS 

We now present examples of the cooperative effects 
of the hydrodynamic interaction between many apolar 
swimmers, for first regular, and then random, distri- 
butions of swimmers. Consider a chain of N identical 
dumb-bells all oriented along the z— direction and ini- 
tially positioned at equally spaced intervals of 5D along 
the y— axis. To optimise the hydrodynamic interactions 
between nearest neighbours the relative phase between 
any two neighbouring dumb-bells is set to 7r/2, increas- 
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FIG. 2: (a) Configuration of a linear chain of dumb-bells and 
a schematic representation of its evolution. The bulk of the 
chain moves faster than the boundary, leading to dumb-bells 
being left behind at the trailing boundary and the formation 
of 'fast pairs' at the leading boundary, (b) The distances 
moved by each dumb-bell in the chain after 100, 000 swimming 
strokes for chains with different numbers of dumb-bells. 



ing in the positive y— direction. Although this is a highly 
artificial configuration, the lack of any rotation greatly 
simplifies the dynamics, allowing for the effect of chang- 
ing N to be clearly quantified. 

We show in Fig. [5] the initial configuration of the chain 
and the distances moved by each of the dumb-bells as 
their total number increases. Two features are partic- 
ularly noteworthy: firstly the behaviour for iV > 2 is 
significantly different from that for N — 2, and secondly 
the general behaviour for TV > 20 shows only minor vari- 
ations, indicating that an asymptotic limit is being ap- 
proached. 

The evolution of the pattern of swimmers may be un- 
derstood by noting that there is a fundamental distinc- 
tion between dumb-bells with two nearest neighbours and 
those with only one. The former constitute what we shall 
call the bulk of the chain, while the latter form the bound- 
ary. A dumb-bell in the bulk not only gets pulled along 
by its neighbour in front of it, but is also pushed by 
its neighbour from behind. These two interactions add 



constructively leading to a member of the bulk moving 
faster than a single isolated pair. By contrast a dumb-bell 
which is on the boundary only has one nearest neighbour 
and hence does not benefit from this added boost. Thus 
the bulk moves faster than the boundary. 

This has the effect of introducing an asymmetry be- 
tween the two boundaries; the dumb-bell on the trail- 
ing boundary gets left behind while that on the leading 
boundary is caught up. From Eq. ([5]) we see that the 
strength of interactions depends on the separation of the 
dumb-bells as r~*, so that as the trailing dumb-bell gets 
left behind, the pull it receives from the bulk rapidly 
diminishes until it becomes isolated and can no longer 
move. At the same time the distance between the lead- 
ing dumb-bell and the bulk decreases yielding a sharp 
increase in the strength of interaction between itself and 
its neighbour. As a result the leading pair speed up sig- 
nificantly and are ejected from the front of the chain. 



LATTICE PUMPS 

The cooperative motion of one dimensional chains car- 
ries over to regular arrays in two dimensions. As an ex- 
ample, consider a square lattice of dumb-bells all oriented 
in the z— direction and with {x,y) positions {jL,kL), 
where L is the lattice constant and j, k are integers. Co- 
operative directed motion can be induced by defining the 
phase of each dumb-bell to be (j)j^k = + k)/2. The 
entire lattice then moves uniformly along the [110] direc- 
tion. 

For a system of real swimmers this state could not 
be sustained; we have found that it is unstable to any 
imperfections of the lattice such as boundaries or fluctu- 
ations in the position of the swimmers. This instability 
of long range coherent states is in agreement with an- 
alytic work by Ramaswamy and co-workers [ill . IT^ . 113 
and numerical simulations of Saintillan and Shelly [l5|. 
However one might envisage tethering fabricated dipo- 
lar swimmers to a substrate and aligning and activating 
them with a magnetic field. Such a set-up would act as 
a micon-scale pump. We estimate v ^ ^ms^^, which is 
similar to the velocities achieved by bacteria. It should 
be noted that the interactions scale with separation as 
r^^, so that a small decrease in lattice spacing will pro- 
vide a substantial increase in flow speed. 



SUSPENSIONS OF DUMB-BELLS 

We consider a suspension of N dumb-bells initially dis- 
persed throughout a cubic box of side L = 201? with 
random positions, orientations and relative phases. To 
avoid singularities in the hydrodynamic interactions a 
short distance cut-off is employed when the separation 
between any two dumb-bells becomes less than 0.51?. At 
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FIG. 3: (Colour online) Collective properties of a suspension 
of interacting (I) dumb-bell swimmers (apolar) and (II) three- 
sphere swimmers (polar). Top; mean speed of the swimmers 
as a function of their number density. The line indicates a 
linear fit. Middle: mean square displacement as a function 
of time for number densities 0.1 (x) and 0.02 (o). Bottom: 
diffusion constant for increasing number of swimmers. 



separations of order D an expansion of the Oseen tensor 
in powers of D /r does not converge rapidly. To overcome 
these difficulties we instead use an expansion in powers 
of ^/r to better describe the near field hydrodynamics. 

The combination of a random initial configuration and 
apolar symmetry means that there is no prefered direc- 
tion for the motion, and on average the velocity is zero. 
In continuum models of apolar active fluids spontaneous 
symmetry breaking can lead to a state with non-zero av- 
erage velocity [l6j . however we have not observed any 
such transitions in our simulations of dumb-bells. In the 
absence of a net velocity, the mean speed provides a mea- 
sure of the degree of collective activity. Fig. shows 
that the mean speed of the dumb-bells increases linearly 
with the number density n = N/L^. This agrees with 
a simple scaling argument: since every dumb-bell will 
interact with every other one, the total number of inter- 
actions scales as n^. Balancing this against the kinetic 
energy predicts that the mean speed should increase lin- 
early with n. The same linear scaling, but tending to a 
finite value as TV — > 1 because of the finite speed of a 
single swimmer, is also found for a suspension of polar 



swimmers (Fig. |3l[II)). 

A measure of the nature of the collective motion gen- 
erated by the interacting dumb-bells is the mean square 
displacement, {Ar'^{t))e = (|r(t) — r(0)p)e, where (•)e de- 
notes an ensemble average. This is plotted in Fig. [3] for 
number densities of n = 0.02 and n — 0.1. At long 
times, and for large number densities, a scaling form 
(Ar^(t))e ^ develops, with an exponent of z = 1 in- 
dicating that the suspension of dumb-bells is behaving 
diffusively. At smaller number densities the behaviour 
is more sporadic because the large average separation 
between dumb-bells greatly reduces the strength of the 
interactions between them. The collective motion is then 
dominated by those fluctuations in the local number den- 
sity which bring two dumb-bells close enough to allow 
them to move appreciably. 

The mean square displacement of a suspension of 
dumb-bells shows qualitatively different behaviour to 
that of a suspension of polar swimmers. In the latter 
case, shown in Fig. [31^11), the mean square displacement 
is ballistic at short times (z = 2), with a cross-over to dif- 
fusion (z = 1) at longer times. Moreover, since this cross- 
over is due to the randomisation of swimmer orientations 
through hydrodynamic interactions, it occurs at later and 
later times as the number density is reduced. For po- 
lar swimmers this leads to a diffusion constant which 
decreases as the number density is increased. For ap- 
olar swimmers the oppposite is true; their motion arises 
solely from hydrodynamic interactions and the diffusion 
constant increases with increasing number density. 

DISCUSSION 

The aim of this letter is to discuss the motion of a 
simple model of apolar swimmers, systems which move 
at zero Reynolds number only in the presence of other 
swimmers. We have demonstrated that, as is the case 
for colloids [13] and polar swimmmers [1] , hydrodynamic 
interactions lead to complex collective behaviour. 

Particular observations are that, for two dumb-bell 
swimmers which are initially parallel, the most likely fi- 
nal state corresponds to the axes of the swimmers ly- 
ing at right angles. Oscillatory trajectories are also ob- 
served. For regular arrays of swimmers cooperation be- 
tween hydrodynamic pair interactions can lead to simple 
flow fields. In particular a square array of dumb-bells 
produces a constant flow and hence, if the dumb-bells 
were fixed in position, could act as a pump. For a large 
number of dumb-bell swimmers, initialised with random 
positions and phases, the mean speed is linear in the num- 
ber of dumb-bells, reflecting the energy pumped into the 
system by the hydrodynamic interactions. The mean- 
square displacement evolves linearly in time showing the 
expected diffusive behaviour. 

There is much further work to the done to explore the 
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phase space of apolar swimmers, for example consider- 
ing initial configurations for which pairs move out of the 
plane, three dumb-bell orbits, and the role of the distribu- 
tion of phases on multi-dumb-bell motion. For more than 
two dumb-bells the relative phase becomes an important 
variable; since the pairwise interactions between three or 
more dumb-bells cannot simultaneously take their maxi- 
mum value the system exhibits a type of frustration. It is 
also interesting to consider the effect of moving away from 
the zero Reynolds number limit , and to ask which of 
the properties of dumb-bells provide a generic representa- 
tion of the class of apolar swimmers. Comparing simple 
microscopic models of polar and apolar swimming may 
help to formulate the correct continuum theory of swim- 
mers, and to link the microscopic and continuum length 
scales. 

We are grateful to Mike Gates, Davide Marenduzzo 
and Chris Pooley for useful discussions and thank Eric 
Lauga and Denis Bartolo for showing us reference 
prior to publication. 
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